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By using renormalization-group (RG) methods, we study a non-mean-field model of a spin glass 
built on a hierarchical lattice, the hierarchical Edwards-Anderson model in a magnetic field. We 
investigate the spin-glass transition in a field by studying the existence of a stable critical RG 
fixed point (FP) with perturbation theory. In the parameter region where the model has a mean- 
field behavior—corresponding to d > 4 for a d-dimensional Ising model -we find a stable FP that 
corresponds to a spin-glass transition in a field. In the non-mean-field parameter region the FP above 
is unstable, and we determined exactly all other FPs: to our knowledge, this is the first time that all 
perturbative FPs for the full set of RG equations of a spin glass in a field have been characterized in 
the non-mean-field region. We find that all potentially stable FPs in the non-mean-field region have 
a nonzero imaginary part: this constitutes, to the best of our knowledge, the first demonstration 
for a spin glass in a field that there is no perturbative FP corresponding to a spin-glass transition 
in the non-mean-field region. Finally, we discuss the possible interpretations of this result, such 
as the absence of a phase transition in a field, or the existence of a transition associated with a 
non-perturbative FP. 
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I. INTRODUCTION 


Establishing the existence of a phase transition in non¬ 
mean-field spin glasses in a magnetic field is an open prob¬ 
lem that has been attracting growing interest in the last 
few decades. Indeed, a demonstration of the existence 
or of the absence of such a transition could shed light 
on a fundamental problem: understanding the structure 
of the low-temperature phase of spin glasses. Namely, 
the occurrence of a transition in a field would indicate 
that the energy landscape of a spin glass has a com¬ 
plex structure involving an exponentially large number 
of states [1], while the absence of this transition would 
hint the existence of only two low-lying energy basins, 
which are reminiscent of the ground states of Ising fer- 
romagnets [2], Despite decades of intense research, the 
occurrence of a phase transition for non-mean-field spin 
glasses in a magnetic field is still under debate: Indeed, 
experimental studies on disordered magnetic materials 
provided evidence both in favor [3] and against [4] a tran¬ 
sition. In addition, numerical simulations for Ising spin 
glasses with short and long-range interactions are affected 
by long equilibration times and strong finite-size effects, 
preventing numerical approaches from giving a definite 
answer on the existence of a transition in a field [5-9]. 

On a conceptual level, one of the reasons why non- 
mean-field models of spin glasses are hard to solve is that 
the complex structure of the low-lying energy states hin¬ 
ders the use of the renormalization-group (RG) coarse- 
graining methods widely employed for homogeneous sys¬ 
tems, such as Ising ferromagnets [10]. In this regard, the 
RG transformation for ferromagnetic spin systems takes 
a strikingly simple form when applied to models built on 
a hierarchical lattice [11]. To study non-mean-field spin 
glasses with RG methods, it is thus natural to consider 
a spin-glass model where spin interactions are disposed 


in a hierarchical way. This model is known as the hier¬ 
archical Edwards-Anderson model [12], and it has been 
raising interest in recent years because it provides a nat¬ 
ural implementation of the RG transformation, allowing 
for a novel RG characterization of the critical properties 
of a non-mean-field spin glass [13-16]. 

In this paper, we study the existence of a phase transi¬ 
tion for the hierarchical Edwards-Anderson model with 
an external magnetic field (HEAM). The hierarchical 
structure of spin-spin interactions for the HEAM implies 
an exact RG equation for the replicated partition func¬ 
tion with fixed overlap which is identical to the RG equa¬ 
tion for the hierarchical Edwards-Anderson model [12]. 
In perturbation theory, this RG equation takes a partic¬ 
ularly involved form: indeed, the presence of the external 
field implies that the replicated partition function is given 
by a combination of twelve overlap monomials [17], result¬ 
ing into a complex perturbative structure. To carry out 
such an involved perturbative expansion, we developed a 
symbolic manipulation tool which allowed us for deriving 
the RG recurrence equations to lowest order in pertur¬ 
bation theory. We then investigated the existence of a 
phase transition in a field by studying the fixed points 
(FPs) of these RG equations. The FP equations are a 
system of polynomial equations of the third degree which 
may have, in principle, a large number of solutions. In 
the parameter region where the HEAM has a mean-field 
behavior we find a stable FP corresponding to the ex¬ 
istence of a spin-glass transition. In the non-mean-field 
region the FP above is unstable, and all other solutions of 
the FP equations need to be studied to establish the exis¬ 
tence of a stable FP. In this regard, to our knowledge the 
full set of perturbative FPs in the non-mean-field region 
of a spin glass in a field has never been fully character¬ 
ized: in particular, two previous RG studies addressed 
this problem for short-range spin glasses in a field. Pi¬ 
mentel et al. found a set of unstable FPs [17], but these 
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FIG. 1. Hierarchical Edwards-Anderson model in a magnetic 
field with k = 3. Dots represents spins Si,... ,Sa from left to 
right, arcs represent interactions between spins below them, 
and magnetic fields are not shown. The red path starts from 
the two circled spins Si, S 3 and it goes from the bottom to 
the top of the hierarchical tree until a common arc between 
the two spins is found: since this requires ascending two hier¬ 
archical levels, the hierarchical distance between Si and S 3 is 
di3 = 2. 

were not shown to coincide with the full set of solutions of 
the FP equations [18]: as a consequence, the possibility 
that other stable FPs could exist remained open. Other 
RG approaches considered only certain linear combina¬ 
tions of the overlap matrix—the ‘replicon’ modes—thus 
retaining only a subset of the twelve monomials in the 
replicated partition function [19, 20]: this approach re¬ 
sulted in a reduced set of FP equations which were found 
to have no stable FP. By using the Grobner-basis method 
for systems of polynomial equations [ 21 ], here we deter¬ 
mined exactly all perturbative FPs of the full set of RG 
equations in the non-mean-field region of the HEAM, and 
we show that all potentially stable FPs have a nonzero 
imaginary part. To our knowledge, this result constitutes 
an unprecedented demonstration that there is no pertur¬ 
bative FP corresponding to a spin-glass transition in the 
non-mean-field region of a spin glass in a field. 

The paper is organized as follows: In Section II we 
introduce the HEAM, in Section IIA we consider the RG 
equations and we discuss their perturbative solution, in 
Section IIB we derive the qualitative structure of the 
critical FP, which is then determined explicitly in Section 
IIC. Finally, Section III is devoted to the discussion and 
interpretation of the results, as well as to an outlook of 
future studies. 


II. RESULTS 


distance [ 22 ], {Jij} are independent and identically dis¬ 
tributed (HD) Gaussian random variables with zero mean 
and unit variance, and the single-spin Hamiltonian on site 
i is H 0 [Si\ = hiSi, where {hi} are IID Gaussian random 
variables with zero mean. Given two spins Si, Sj, let us 
consider the number of levels that we need to ascend in 
the hierarchical tree starting from the bottom in order to 
find a root common to Si and Sj, see Fig. 1: we denote 
this number by the hierarchical distance dij. The recur¬ 
sive structure of the Hamiltonian (1) then implies that 
the interaction strength between spins Si, Sj depends 
on i, j through their hierarchical distance rather than 
through their Euclidean distance: indeed, Eq. (1) can 

be rewritten as H k [S\ = - Y?i<j=i K ijSiSj + Yh=i h i S i> 
where {Kij} are IID Gaussian random variables, and K l3 
has zero mean and variance ( YaZq ' 3 2 ~ 2 al ) 2 ~' 2 ' jdij . 

To relate the HEAM to other spin-glass models with 
different disorder distributions such as short- and long- 
range Ising spin glasses [23, 24], two considerations are 
in order. First, a precise correspondence between short- 
range spin models in d dimensions and one-dimensional 
long-range models is still subject of research [ 8 , 25]. Sec¬ 
ond, the universality of physical observables with respect 
to the form of the bond distribution in Ising spin glasses 
is also a subject of ongoing debate [26-29]. Despite the 
considerations above, a hierarchical Edwards-Anderson 
model with a given value of er has some features in com¬ 
mon with a short-range Ising spin glass on a hypercubic 
d-dimensional lattice [13]. In particular, the parameter 
region 1/2 < a < 2/3, where the hierarchical Edwards- 
Anderson model has a mean-field behavior, is analogous 
to d > 6 for a short-range Ising spin glass, i.e. it is charac¬ 
terized by vanishing order-parameter fluctuations. Simi¬ 
larly, the non-mean-field region 2/3 < a < 1 corresponds 
to d < 6 for an Ising spin glass, and it is characterized by 
a non-Gaussian RG FP associated with nonzero order- 
parameter fluctuations. 


A. Renormalization-group equations 


The hierarchical structure of the interactions in the 
HEAM shown in Fig. 1 implies an exact RG equation 
which can be derived with the replica method [13]: we 
introduce an integer n, an n x n matrix Q a b with Q a b = 
Qba and Q aa = 0 , and the replicated partition function 
with fixed overlap 


The HEAM is a system of 2 k Ising spins S = 
Si,... ,S 2 k, Si = ±1, whose Hamiltonian H k [S] is given 
by the recursive equation [ 12 ] 


Z k [Q}= E 


£e-0E:.i **[«“] 


a<.b— 1 


( 2 ) 


H k [g\=H 1 k _ 1 [Si]+H 2 k _ 1 [S 2 £ JijSiSj, ( 1 ) 

i<j =1 

where Si = Si ,..., S 2 k-i and S 2 = S 2 k-i + 1 ,..., S 2 k de¬ 
note the spins in the left and right half of the lattice 
respectively, 1/2 < a < 1 is a parameter which de¬ 
termines how fast spin-spin interactions decrease with 


where in Eq. (2) E[] denotes the expectation with respect 
to all random variables, /3 = 1/T is the inverse temper¬ 
ature, {S a } are n replicas of the spin configuration S, 
and 

g aiS l£S‘St 


P) 
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is the overlap [1] between replicas a and b. As we will 
show in what follows, the replicated partition function 
(2) is a useful quantity because, in principle, its n —» 0 
limit incorporates all thermodynamical properties of the 
system [1]. Let us introduce 

C = 2 2(1 - ct) (4) 


and the rescaled overlap distribution 

Z k [Q\ = Z k [C- k ' 2 Q]. (5) 

By using Eqs. (1), (2), it can be shown [12, 13] that Z k 
satisfies the following recursive equation 


Z k+ i[Q\ = exp 




Q + P 
C V 2 


Zk 


Q-P' 

CM 2 ’ 


( 6 ) 



TABLE I. Monomials contributing to the overlap probability 
distribution of the hierarchical Edwards-Anderson model in a 
magnetic field to order Q A . The monomials of order one, two 
and three are listed in the rows of the first, second and third 
column respectively, from top to bottom in increasing order 
of the index p. 


where J[dP] = J Y[ a<b dP ab , and in what follows we will 
omit Q-independent proportionality factors multiplying 

Zk- 

Being a functional recursive equation, Eq. (6) is hard 
to solve with exact methods: still, an analytical solu¬ 
tion can be worked out in perturbation theory. We write 


Z k [Q] as a combination of monomials in Q , where each 
monomial is multiplied by a numerical coefficient which 
depends on k: then, Eq. (6) implies a set of algebraic re¬ 
cursive equations for these coefficients, and we will show 
that these equations can be studied analytically in per¬ 
turbation theory. Neglecting monomials of order larger 
than Q 3 , the most general form of Z k [Q] is given by [30] 


Z k [Q] = 


exp 


1 3 
1 \ ' 


- s k l 1 [Q} + -Yrlll[Q} + - 


2 ^ 

p=i 


1 8 

-y , 

6 ^ 

p= i 




u k p 


( 7 ) 


where the monomials I 1 [Q], {Ip[Q]}, {Ip[Q]} are listed 
in Table I. In Eq. (7), Z k is given by a combination of 
twelve terms: given that these are all the monomials 
of degree three or less that are consistent with the 
symmetries of the model (1), Eq. (7) constitutes the 
most general form of the third-order expansion of Z k 
[30]. In the limit of zero magnetic field additional 
symmetries appear, and the set of admissible monomials 
is reduced [24] to l\[Q\ and I 3 [Q]. 


Plugging Eq. (7) into the RG equation (6), we obtain 


Z k + i [Q\ = exp 



x 



2 s k 
CM 2 


I'M 


2^ C 

p =1 


2 ^ ^ Pib' y ^'ab,cd[Q\Pzd j 5 

a<.b,c<Ld J 



( 8 ) 


where ^ a 6,cd is a n(n — l)/2 x n(n — l)/2 matrix given by 

*^'ab,cd[Q\ ^ab,cd “I” ( ^ab,cd\Q\’> 

and 


4 r} 

ab,cd = (6 a 


•6bd + 6b c 6ad) H- (6a 


a r 3 

+ 6bc + 6ad + 6bd) H- 


*3ab,cd[Q\ = 


1 12 wM 

0 | ^3^2 MCid'eJbc ' Obd'°Ja 


+ 5 

ac Qbd 4“ $bcQad ) 


24 u ’ k (§ 
(J3/2 1 a 


z6bd A 6 bc 6 ad ) Q cd 


( 9 ) 


( 10 ) 


( 11 ) 
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+ 


Awi 

C 3/2 


o 4 

(J3/2 [(dacfibd + &bc&ad) I* [Q] + 2 (Q a b + Q cd)] 


+ 


Awl 


{S ac Sbd T Sbc&ad) ^ ' (Qce "F Q de) T (S ac 4“ Sb c 4“ Sad 4” Sbd) Qcd 4” 

e—1 

+S 

ac Qbc &bcQac “I” &adQbd “1“ &bdQad 

n 

5> ac Qde ^bcQde H - ^adQce ^bdQce “I - ^ acQbe H - &bcQae “I” ^adQbe “1“ &bdQae) 

12 w 6 x n A 

/2 ^ acQce “I - &bcQce “I - &adQde “I” &bdQ de) H” 
e=l 

n 

(<5ac + Sbc + Sad + Sbd) I 1 [Q] + 2 (Qae + Qbe + Qce + Qde) 


e—1 


c 3/2 

+Q ac H” Qbc H” Qad H” Qbd 


+ 


c 3 / 2 


e—1 




Note that in Eq. (9) we split the matrix .jW, into a 
Q-independent term & and a Q-dependent term Sf, the 
latter being proportional to the coefficients { w £}: indeed, 
in what follows we will make the assumption that the co¬ 
efficients {ref} in Eq. (7) can be treated in perturbation 
theory, and the form (9) will be convenient for solving the 
RG equation ( 8 ) with an expansion in powers of {u^}. To 
set up this expansion, let us first compute the inverse of 
&. This can be sought as a combination of the terms 
that appear in Eq. (10): 


«•-1 -- 
^ ab,cd 

Cl (fiac^bd H - $bc$ad) ^2 ($ac “1“ $bc $ad &bd) H - ^3,(12) 


and the coefficients Ci, C 2 , C 3 can be determined by solving 
the equation Y, c <d^ab,cd&tiL, e f = S ae Sbf■ The result is 


Cl = 


C_ 

M-’ 


(13) 


__ 

4r fc ( 2r l + (n- 2 )rl) ’ 

C(( r ic) 2 “ 2 r l r k + n r t r k) 

2 r\ ( 2 r 4 +(n- 2 )r 2 ) ( r\+{n - 1 ) (r^+nr 3 ))' 


(14) 

(15) 


The matrix & 1 can now be used to compute the second 
factor in Eq. ( 8 ): 


/ \(1P] exp ( - ^ ^ Pab^ab,cd[Q] P cd\ = 

\ a<b,c<d / 

a 

y/det(^[Q]) 

a exp ^Trlog (l + • Sf[<2])^ = 


a exp I — — 4 V 


^- 1 - S^[Q] - - (^- 1 -^[Q]) 2 + 


+ -(^- 1 -5f[0]) 3 + O(Q 4 ) 


(16) 


where in Eq. (16) a denotes a Q-independent factor, and 
in the third line the symbol ‘ 1 ’ in the logarithm denotes 
the n{n — l )/2 x n(n — l )/2 identity matrix. 

Given the complex form (11) of the matrix Sf, the cal¬ 
culation of the matrix products and traces in Eq. (16) 
is involved: for example, if we denote by a ‘term’ in 
Eqs. (10), (11) the sum of multiple addends related 
to each other by the permutation a 6 , c <-> d — 
such as S ac + Sbc + Sad + Sbd —then a direct evaluation 
of Tr[(^ -1 ^) 3 ] involves the computation of (3 x 14) 3 = 
74,088 matrix traces. To handle this computation, we 
developed a symbolic manipulation program [31] that en¬ 
abled us to compute explicitly the right-hand side (RHS) 
of Eq. (16): Sums over replica indexes are computed in 
a symbolic way, and the result depends explicitly on the 
number of replicas n, allowing for a straightforward ana¬ 
lytic continuation [1] for n —> 0. Importantly, the result 
of this calculation shows that the RHS of Eq. (16) is 
a linear combination of the monomials listed in Table I. 
This implies that -Zk+ilQ] and ^k[Q\ are both given by 
a linear combination of the same monomials, hence the 
RG equation ( 6 ) can be solved consistently in perturba¬ 
tion theory. In particular, Eqs. (7), ( 8 ), (16) imply a set 
of recursive equations which relate Sfc+i, {r^ +1 }, {w^ +1 } 
to Sfc, {r^}, {w^}: the n —> 0 limit of these equations is 
given by Eqs. (S1)-(S12) in Section SI of the Supplemen¬ 
tal Material. In addition, in Section S2 we discuss how 
to write the RG coefficients s*,, {r^}, in terms of 

the interaction-decay exponent, the temperature and the 
magnetic-field strength to study, for example, how the 
RG flow or the FP domains of attraction are related to 
the physical parameters of the model. 


B. Structure of the critical fixed point 


We will now study the existence of a FP associated with 
a spin-glass transition and determine its general features. 
To do so, we consider the spin-glass susceptibility [ 8 , 9] 

2 ^ 

XSG=^J2 - (SiHSj)) 2 }, (17) 

i,j =1 
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where () denotes the Boltzmann average with Hamilto¬ 
nian (1). The susceptibility (17) can be rewritten as [32] 

Xsg = 2 fc E[(g? 2 ) _ 2(912913) + (912934)], ( 18 ) 

where in what follows the average () of a function of 
multiple replicas denotes the average over all replicas [1], 
each replica having an independent Boltzmann measure 
with Hamiltonian (1). By using Eq. (S20), the overlap 
averages in the RHS of Eq. (18) can be rewritten as 
follows: 


E[<<&)] = <?"* lirn Qf 2 , 

n—> 0 

(19) 

£[(912913)] = C k lim Q 12 Q 13 , 

n—>-0 

(20) 

£[( 912934 )] = C k lim Q 12 Q 34 , 

n—>-0 

(21) 


where denotes the average with respect to Z k [Q\. Ne¬ 
glecting 0(w) terms and using Eqs. (7), (10), we obtain 


Zk[Q] = 


exp 


C 


$k I [Q] T , ^ ( Qab^ab,cdQcd 


a<.b,c<d 


• ( 22 ) 


By using Eq. (22) and standard Gaussian-integration 
rules [33], we have 


Qi 2 


Q 12 Q 13 


Q12Q34 


(^) ^ 2 + ^r 2 p 2 ( 23 ) 

(~cr) + c^ Cl + 2c2 + C3 )’ 

(^) ^ 2 + §^r 2 :i 3 (24) 

+ (o( C 2 + c 3)> 

(^f) ( 25 ) 


would imply that xsg is finite. We recall that the con¬ 
clusion above on the large-fc behavior of r\ at the critical 
point has been derived by neglecting the 0(w) terms in 
Eq. (7): given that our RG analysis is based on the work¬ 
ing hypothesis that perturbation theory is well-behaved, 
retaining the 0(w) contributions would result in a pertur¬ 
bative correction to r), without changing the qualitative 
behavior (27). 

We have thus derived a first property of the critical 
FP: at the critical temperature T = T c , we have r l —» r\ 
for k —» 00 , with r* finite. Importantly, Eq. (S2) shows 
that if T ^ T c , then r/ diverges like r\ ~ (2 /C) k for large 
k: setting T = T c , the projection of the RG flow along 
the direction associated with r\ is set to zero, and the 
divergence of r\ is removed [10]. 

The critical FP can now be completely characterized 
as follows. According to Eqs. (SI), (S3), (S4), the 
coefficients Sk,rl,r\ diverge for large k: we cannot re¬ 
move these divergences like we did for the coefficient 
r\, because the only free parameter in the model -the 
temperature—has already been fixed to eliminate the di¬ 
vergence of rl. Thus, the only possible FP for these coef¬ 
ficients is Sfc = s* = ± 00 , = rl = ± 00 , r\ = rl = ± 00 . 

Finally, given our perturbative working hypothesis, the 
coefficients {u^} are assumed to be small, thus they must 
converge to a finite FP = u>*, p = 1,..., 8. 

Given the critical-FP structure above, the critical val¬ 
ues {c*} of the coefficients {c p } are finite despite the fact 
that s*, rl and rl are infinite: 


C 

4rp 

C 

8^’ 

C 


4 rl 


(28) 

(29) 

(30) 


where to obtain Eqs. (28)-(30) we used Eqs. (13)-(15) 
and the conditions rl = ± 00 , rl = ± 00 , n —> 0. Impor¬ 
tantly, the finiteness of {c*} ensures that the FP equa¬ 
tions for rl, {ref} are well posed because they involve 
only finite terms—see Eqs. (S2), (S5)-(S12). 


where in Eqs. (23)-(25) we set Jif = Yh a <b ab an< l 
we used Eq. (12). Putting together Eqs. (4), (13), (18), 
(19)-(21), (23)-(25), we obtain the expression for the spin- 
glass susceptibility as a function of r\\ 

( 2 \ / 2 \ Q 

c) tv i = {c) (26 > 

Since we are assuming that xsg diverges at the critical 
point, if any, the coefficient r\. must have the following 
behavior 

r\ r\ for k —» 00 , T = T C with r\ finite. (27) 

Indeed, if Eq. (27) did not hold, then according to Eq. 
(S2) rj. would diverge like rl ~ (2 /C) k , and Eq. (26) 


C. Solution of the fixed-point equations 


We will now determine explicitly the critical FP. First 
of all, the RG Eqs. (S2), (S5)-(S12) possess a trivial FP 


P 2 

- W 

2{2/C — 1) ’ * 


= wl = 0. (31) 


As shown in Appendix A, in the mean-field region 1/2 < 
a < 2/3 this FP is stable, and it is thus associated with 
the existence of a physical spin-glass transition in a mag¬ 
netic field. 


In non-mean-field region 2/3 < a < 1 the FP (31) is 
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unstable, see Appendix A, thus other FPs must be consid¬ 
ered. To this end, we observe that the FP equations can 
be simplified if er lies in the neighborhood of the threshold 
value cr = 2/3. To show this, we set 

e = er - 2/3, (32) 


(33) and the recursive equation (S2) imply that the FP 
for r* reads 


r 


l 

* 


( 3 2 

2(2V3_l) +pe ’ 


(34) 


and we observe [13] that for small e the RG Eqs. (S2), 
(S5)-(S12) are of the form w*e = 0((w *) 3 ), implying that 
( w S) 2 = O(e). Hence, we set 


„P - P 3 


log 2 

(2l/3_l)3j 


1/2 


(33) 


where {w p } are coefficients of order unity [34]. Also, Eq. 


where p is a coefficient of order unity. Note that the 
mean-field FP (31) can be obtained from Eqs. (33), (34) 
by setting e = 0: hence, the FP (33), (34) can be re¬ 
garded as a perturbation of the mean-field FP. We will 
now determine the critical FP, if any. To this end, we use 
Eqs. (33), (34), (28)-(30), and we obtain [31] that the FP 
condition for Eqs. (S2), (S5)-(S12) implies the following 
set of polynomial equations for p, {ui p } 


4(2 1/3 - l) 2 p + 2 1/3 /3 2 log 2 (-4 uj\ A 16wiw 2 - Wu\ + 4) = 0, (35) 

14w 3 — 36w 2 ca 2 T 18cuica 2 -f- T 6uq = 0, (36) 

w 2 (3 + 12 uf - 30wiw 2 A 17 ul) = 0, (37) 

18w? — 4w 2 (15w 2 + 2^3) A 8wiw 2 (9w 2 + 4^3) — llw|(3w 2 + 2w 3 ) — 6W3 = 0, (38) 

6 uj 3 A 8ca 2 (c^3 A 2uq) — aqca 2 (15ca 2 A 12c^3 A 64uq) A uj 2 (12cu 2 A 7cu 8 A 44uq) A 12aq = 0, (39) 

18w 3 — 3wj(27w 2 A 4ca3) A 6uq (6c^ 2 — 3uj 2 ca3 — cu 8 ) A ca 2 (18(^2 A 24cu 2 ui3 A Tcu^'j A I8UI5 = 0, (40) 

72uj 3 — 243 uJ 1 uj2 A 18cui (l2cu2 A ca 2 LU 3 A cc 8 ] — 63ce 3 — 18ca 2 ca3 — 15 w 2 cu 3 A 4ca 8 — 54ca 8 = 0, (41) 

9w 3 A 3uq (6 w 2 A 7w 3 A 8uq) A 3uq [3ca 2 A 4cu 2 (ca3 A 4uq) A cu3(3ca3 A 16uq)j — 36w 2 — 3(^2(llu?3A (42) 

A 2 OW 4 ) — 2w 2 W3(5ciJ 3 A 24uq) — w 2 (w 3 — 8 W 4 ) — 36W7 = 0, 

63w 3 — 9 w 2 [3w 2 — 5(w 3 A 4 ^ 4 )] — 9wi(w3 A 4 ^ 4 )(2w 2 — W 3 — 4uq) — 36w| — 36wf (W 3 A 2uq) — 3w 2 x (43) 

x (5w 2 A 24w 3 w 4 A 48 (^ 4 ] — 0 J 3 A 12 cc 3 uq — 48w 2 ) — 216 uj 8 = 0. 


The solution to the equations above can be determined 
by solving Eqs. (36)-(43) for {lj p } first, and then substi¬ 
tuting the solution into Eq. (35) to obtain p. Being a 
system of polynomial equations of the third degree, in 
general Eqs. (36)-(43) have multiple solutions. To find 


all solutions, we determined the Grobner basis (GB) of 
the system of polynomials (36)-(43): the GB is given by 
a set of polynomials in {w p }, and the roots of Eqs. (36)- 
(43) coincide with those of the GB [21]. Given that Eqs. 
(36)-(43) have integer coefficients, their GB can be com¬ 
puted exactly, and it reads 


w 8 (28w| A 3) {16 [80 (8656w| A 933) uj A 2367] vj A 3267} , (44) 

w 7 - 3w 8 , (45) 

1877310171w 6 A 32100226170880wJ A 2352496281600w| A 641548992w| A 8790294330w 8 , (46) 

625770057w 5 - 8w 8 (28wf A 3) (71652290560*4 - 2425923360w| A 261352389), (47) 

625770057^4 A 8025056542720wJ A 588124070400w| A 160387248w| A 1258918497w s , (48) 

625770057w 3 A 2w 8 (8025056542720w| A 588124070400*4 A 160387248*4 A 1258918497), (49) 

1877310171w 2 - 16w 8 (28*4 A 3) (71652290560*4 - 2425923360*4 A 261352389), (50) 

625770057cji - 2w 8 (8025056542720*4 A 588124070400^ A 160387248*4 A 2510458611), (51) 


where every line in Eqs. (44)-(51) corresponds to an element of the GB. The first GB element, Eq. (44), 
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depends only on w 8 , and its complete set of roots 
can be determined exactly. Since the GBs (45)-(51) 
are linear in lo \to every solution for w 8 cor¬ 
responds a unique value for uq,..., W 7 . As a conse¬ 
quence, the GB method allows for extracting the full 
set of exact solutions. The numerical root values for 
w 8 read w 8 = 0, 0.159676 ± 0.167743 i, -0.159676 ± 
0.167743 i, ±0.327327 i, ±0.320162 i. First, it is straight¬ 
forward to show that the trivial root w 8 = 0 corresponds 
to the FP 

„ 2 1 ' /3 1os'2 

P = ~P ( 91/3 _ 1 ^ 2 ’ ^i = ---=^8 = 0, (52) 

and that this FP is unstable, thus it does not correspond 
to a spin-glass transition—see Appendix B for details. 
Second, all other roots for w 8 have a nonzero imaginary 
part. This implies that, since the initial condition Z 0 [Q\ 
of the RG transformation ( 6 ) is real and since the RG 
transformation maintains reality, none of the FPs in the 
non-mean-field region a = 2/3±e are physically accessible 
by the RG flow. 


III. CONCLUSIONS 


Establishing the existence of a phase transition in non¬ 
mean-field spin glasses with an external magnetic field is 
an open problem which has been attracting growing inter¬ 
est in recent years [3-9]. Indeed, the occurrence of such a 
transition is believed to be related to the structure of the 
low-temperature phase of spin glasses—a central topic 
in statistical physics of disordered systems [2]. Among 
non-mean-field models of spin glasses, the hierarchical 
Edwards-Anderson model [12] is the simplest non-mean- 
field spin-glass system where the hierarchical structure 
[ 11 ] of spin interactions allows for a natural implemen¬ 
tation of renormalization-group (RG) techniques: recent 
studies showed that these RG methods provide a novel 
way of understanding the thermodynamical properties of 
the model [13, 15, 16, 35]. 

In this paper, we studied the existence of a phase tran¬ 
sition for the hierarchical Edwards-Anderson model with 
an external magnetic field (HEAM): We used a RG ap¬ 
proach based on the replica method [13], and we ana¬ 
lyzed in perturbation theory the RG flow of the average 
replicated partition function with fixed overlap. In this 
approach, the replicated partition function is described 
by a set of twelve parameters: given the complex alge¬ 
braic structure of the resulting RG theory, we developed 
a novel symbolic computation method which allowed us 
for extracting the RG equations to lowest order in per¬ 
turbation theory. 

In the absence of a magnetic field, the HEAM pos¬ 
sesses a mean-field region 1/2 < a < 2/3 where order- 
parameter fluctuations vanish, and a non-mean-field re¬ 
gion 2/3 < cr < 1 characterized by nonzero order- 
parameter fluctuations, where cr is a parameter tuning 
the interaction decay with distance [13]. We investigated 
the occurrence of a spin-glass transition in the HEAM by 


studying the existence of a stable RG fixed point (FP) 
in both these parameter regimes. In the mean-field re¬ 
gion we found a stable critical FP which corresponds to 
a physical spin-glass transition. In the non-mean-field 
region, we investigated the FPs contiguous to the mean- 
field one with a perturbative expansion in e = cr — 2/3. 
The resulting FP equations are a system of polynomial 
equations of the third degree which possess, in princi¬ 
ple, multiple FP solutions. In this regard, two previous 
RG studies investigated the perturbative FPs in the non¬ 
mean-field region of short-range spin-glass models in a 
magnetic field. First, a set of unstable solutions of the 
FP equations was found [17], but these solutions were not 
shown to be the complete set of roots of the FP equa¬ 
tions [18]. Second, the complete set of FPs has been 
extracted in a reduced framework where only a subset 
of the order-parameter modes is retained [19, 20]. By 
using a Grobner-basis method for systems of polynomial 
equations [ 21 ], here we computed exactly all perturbative 
FPs for the full set of RG equations in the non-mean- 
field region of the HEAM. Our analysis shows that that 
all potentially stable FPs have a nonzero imaginary part. 
Given that in the RG transformation the initial values 
of the coefficients are real, and given that the RG trans¬ 
formation maintains reality, our results constitute, to the 
best of our knowledge, the first demonstration for a spin 
glass in a field that the there exists no perturbative FP 
in the non-mean-field region. 

To interpret the absence of a perturbative FP in the 
non-mean-field region, several scenarios can be considered 
[19]. First, the absence of a perturbative FP may imply 
that there is no spin-glass transition in the non-mean- 
field region [ 8 , 20]. In this regard, we recall that the 
addition of a random magnetic field in the ferromagnetic 
Ising model increases the lower critical dimension from 
d = 1 in the zero-field case to d = 2 in a finite field [36]: 
along these lines, the inclusion of a random field in the 
hierarchical Edwards-Anderson model may decrease the 
value of a corresponding [ 8 ] to the lower critical dimen¬ 
sion from cr = 1 in zero field to cr = 2/3 in a finite field, 
thus providing a possible explanation for the absence of 
a transition for 2/3<cr<l. A second possibility is that 
a FP in a field for cr = 2/3 ± e exists, but this FP may 
not be contiguous to the mean-field one, thus it may not 
be found with the perturbative approach used here: this 
hypothesis is in line with the fact that perturbative ap¬ 
proaches are generally not well understood in spin glasses 
[13, 15, 37]. A third possibility is that there is a phase 
transition in a field, but this transition is not associated 
with a FP within the replica RG method. As a future 
direction, the last two possibilities may be investigated 
with recent real-space RG approaches developed for the 
hierarchical Edwards-Anderson model [15, 16] which do 
not rely on perturbation theory, nor they make use of 
the replica formalism. Non-perturbative effects could be 
also studied with Monte Carlo (MC) simulations. In this 
regard, recent MC studies focused on four-dimensional 
short-range spin glasses [7] and one-dimensional long- 
range spin glasses which correspond to short-range mod¬ 
els with d = 4 [9]. These works hinted at the existence of 
a transition in a magnetic field in the non-mean-field re¬ 
gion, a picture which is at variance with the perturbative 



results provided here for the HEAM. Non-perturbative 
effects could be pinned down directly in the HEAM by 
means of a systematic comparison between the pertur¬ 
bative RG predictions provided here and MC simula¬ 
tions [38]: indeed, the RG flow of the replicated parti¬ 
tion function Zk —> Zk+i could be directly investigated 
numerically by computing the moments of the overlap 
distribution for different system sizes 2 k . If the numer¬ 
ics provided evidence for a transition in the non-mean- 
field region, one could then probe directly the nature of 
the FP corresponding to such transition by characteriz¬ 
ing numerically its critical exponent v [10]. If the MC 
estimate of v was found to be close to the classical value 
= 1/ (2a — 1) [13] for a > 2/3, then the spin-glass 
transition resulting from the numerics could be associ¬ 
ated with a FP which is a perturbation of the mean-field 
one. Conversely, a strong discrepancy between the MC 
estimate of v and u c \ for a > 2/3 would hint at the exis¬ 
tence of a non-perturbative FP lying outside the domain 
of attraction of the mean-field FP. Finally, a further nat¬ 
ural way to examine non-perturbative effects consists in 
studying the behavior of the e-expansion to large orders 
[13]: to this end, the symbolic computation method in¬ 
troduced in this paper may be helpful in setting up a 
fully automated, large-order e-expansion, which may be 
useful for understanding the limits of perturbation theory 
in spin glasses. 


transformation Xk —t Xk +i implied by Eqs. (S2), (S5)- 
(S12) in the neighborhood of the FP cc* = (r*, 0,..., 0): 
the FP is stable if the matrix 


A4ij — 


dxi 


k +1 


dxl 


(Al) 


has not more than one eigenvalue larger than one [10]. 
By using Eqs. (S2), (S5)-(S12), (Al), it is straight¬ 
forward to obtain the eigenvalues of A4, which read 
Ai = 2/C, A2 = ••• = Ag = 2/C 3 / 2 . It follows that in 
the mean-field region 1/2 < a < 2/3 only Ai is larger 
than one, and the trivial FP (31) is stable. In the 
non-mean-field region A2, • • • , Ag are all larger than one, 
and the FP (31) is unstable. 


Appendix B: Instability of the trivial fixed point (52) 


We will show that the FP (52) is unstable by proceed¬ 
ing along the lines of Appendix A: since we want to study 
the RG flow in the neighborhood of a FP of the form (33), 
(34), we set 
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P 2 


w k = P 


2(2V 3 - 1) 

log 2 


PfcC 

-,1/2 


( 2 1 / 3 - 1)3 


^p kV~e. 


(Bl) 

(B2) 


We then introduce the vector y/ = {y\,...,y\) = 
(pfc, w 1 fc,..., u>8 k) and the RG transformation y/ —> yk+i 
implied by Eqs. (S2), (S5)-(S12), (Bl), (B2). Then, we 
consider the FP y* = (p,u> 1, • • • , w§), where p,u> 1, - ■ • , 
are given by Eq. (52): this FP is stable if the matrix 




dyj+i 

dy{ 


(B3) 


Appendix A: Stability of the trivial fixed point (31) 
in the mean-field region 


To study the stability of the FP (31), we set Xk = 
(xl ,..., x%) = ( Tk , Wf., • • • , tu 3 ) and we linearize the 


has not more than one eigenvalue larger than one. By 
using Eqs. (S2), (S5)-(S12), (Bl), (B2), (B3) we obtain 
the eigenvalues of A f, which read Ai = 2 1 / 3 , A 2 = • • • = 
Ag = 1 + 3 e log 2: given that in the non-mean-field region 
e > 0, the FP (52) is unstable. 
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SI. RENORMALIZATION-GROUP EQUATIONS 


To obtain the recursive equations relating st+i, {r|! +1 }, {u>^ +1 } to s k , { r k }, we proceed as follows. We 

consider Eq. (8), we replace the left-hand side (LHS) with the expression for Zk +i given by Eq. (7), and we substitute 
Eq. (16) into the RHS: In the resulting equation, we match the coefficients of the monomials I 1 , {1^}, {I 3 } in the LHS 
and RHS, and we obtain a set of recursive equations which relate the coefficients at the k + 1-th step to the coefficients 
at the fc-th step. This calculation has been performed with a symbolic computation tool, and it is detailed in the 
online Mathematica [1] notebook symbolic_computation.nb, which is available as an ancillary file. The resulting 
equations depend analytically on the number of replicas n, thus we could easily extract their n —i 0 limit, which reads 
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+24 w 3 k w\wl + 24w 3 k wlw% + 24w^ (^)“ + 24^w^^ - 8 (wl) 3 - 12 (wlf wl] + c? [c 2 ( - 216 (w^) 3 
+9 («4) 2 (3w^ — 42w| + 64w^ + 88w^ + 60w® — 32 w k ) + 6wfc(l8 (w^f — 3 wl(w k + 8(6 w k + 5 wl + 3 wl 
~4 wl)) — 2(25 (wl) — 2wl(26wl + 34 w k + 21w® — 8 w k ) + 6(4 w\(w\ + w®) + wl(5wl + 6u>®)))) 

+18 (u^.) 2 (10w® + 16u^ + 8 wl + 3w k — 16w^) + 3wl(37 (w k f — 8wl(22wl + 19 w k + 12w® — 8 w k ) 

+12(8 wf,(wl + wl) + wl(5wl + 6u>®))) — 2(33 (w|) — 2 (w^) (36+ 46w® + 27w® — 8 w k ) 

+24wl(2wi + 3 w 5 k )(wl + - 12 (w 5 k f (2wl + 3^))) + 2c 3 (3wl - 3w 2 k + w 3 k ) 2 (&w l k + 7wf. - 6(w^ 

+«>fc))] + 2 cic 2 [c 2 (486 (wfc) 3 - 9 (wlf (57- 96u| + 48w^ + 104^ + 78w| - 16 wl) + 6w^(9 (u^.) 2 
+3wl(-45wl + 40w^ + 56w®. + 42^ - 16w|) + 73 (w 3 k f - 2w 3 k (28wt + 70w| + 45w^ - 8 w 7 k ) + 12(2^ 
x ( wl + w t) + w U 5w k + 6 w 6 k ))) + 54 (wlf + 18 (w 2 k f (9wl - 4(4 w^ + 4 w 5 k + 3 wl - 2w 7 k )) - 3w 2 k 
x (79 (wlf - 4w 3 k (24wi + 40m^ + 27 wl - 8 w 7 k ) + 12(4 wt(w\ + w%) + Wfe(5w| + 6^))) + 2(34 (wlf 
- ( w lf ( 32w fe + 88w l + 51w fc - 8w l) + 24 w 3 k (wl + 3^)(w| + tuf) - 12 (w 5 k f (2w 5 k + 3u^))) - 12c 3 
x (3 wl - 3 wl + wl) 2 (2wl — wl+wl - — to®)] + 4c2(6w^. — 3wl + 2to^ — 2wl) 2 (—6wl + 3tOfc 

—4^ 3 - + 4toj( + 6to®)}, 

Oo/;4 0/1 o 9 9 

+ 2 7(g9/2 {[ ( w fc) + 2 (10to^ + 3(to| + to® — 3to^ + 4to®)) (wl) + 2( — 36 (to^) - 4(7w\ (S8) 

+3to® — 4wl)wl + wl(8w k — 7w\) + 6wl(4wl + w k + 3toJ - 12to®))to^ + 9 (wl) 2 (wl - 4w% + 4to^ 

-2to| + 24t ul) + 4(9(to^ + to| + 6to|) (to^.) 2 + 3(l4 (u^) 2 + 6to|to^ + (tof.) 2 - 2(4wi + wl)w 7 k )w 2 k 
+2(3 (w 5 k + to£) (4) 2 + (wlf wl + (to]!) 2 (to| - w 7 k ))) + 3«4( (to],) 2 + 4(t0] - 5to] + 12to])to] 

-2(32 (to]) 2 + 16to]to] + 7 (to]) 2 - 8(2to] + w 5 k )w 7 k + 6to](to] + 12to])))]c? + [c 2 (81 (wlf + 9(5to] 

—4(7to] + 8to] — 14to] + 24to])) (wlf + 3(9 (to]) 2 — 8(23to] + 12to] + 3to] — 22to] + 24to])to] 

+4(56 (to]) 2 + 52to]to] + 12to®to] + 23 (to]) 2 — 8(4to] + 3to])to]) + 12to](to] — 2(to] + 8to] 

—24to])))to] + 15 (to]) 3 - 48 (to]) 3 + 256to] (to]) 2 + 92to] (to]) 2 - 48to] (to]) 2 - 188 (to]) 2 to] 

—88 (to]) 2 to] — 96 (to]) 2 to] + 304to]to]to] — 48 (to]) 2 to® — 96 (to]) 2 to] + 144to]to]to] + 120 (to])” to] 

+48 (to]) 2 to] — 128to]to]to] — 96to]to]to] + 36 (to]) 2 (6to] + 2to] + 5to] — 24to]) — 96 (wl)" wl 
+12to]( (to]) 2 + (22to] + 8to] + 6to] — 32to] + 48to])to] — 4(12 (to]) 2 + (10to] + 3to® — 8to])to] 
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+ w k( 2w k - 3 w 7 k )))) - 4c 3 (3«4 - 3 w 2 k + w 3 k )(2 (wl) 2 - {llw$. + 5+ 3 w% - 4w 7 k )w 3 k + 6 wl(w 5 k + w%) 
+3wl(7wl + 2w\ — 4w 7 k ) + 3wl(wl - 9 w\ — 2w\ + 4w 7 k ))]cl — 2[(135 (wl) 3 + 9(23 wl — 4(15 
+15w k + 3 w k — 17 wl + 12io®)) (wl) 2 — 3(18 (wl) 2 + 12(2io® — 14 w k — 11 w k — 6 to® + 22w 7 k — 24 w k )w k 
-47 (wl) 2 - 4(24 (wl) 2 + 44 w\w\ + 12 'w k w k + 25 (wl)" — 8(2 wl + 3w\)w 7 k ) + 8wl(25wl + 20 wl 
+ 6 w% - 21 w 7 k + 12 wl))w l k + 33 (wlf - 48 (wl) 3 + 112to| (wl) 2 + 100to| (w 5 k ) 2 - 48 wl (w 5 k ) 2 

— 156 (wl)" wl — 112 (wl) 2 wl — 48 (wl) 2 wl + 224 wlwlwl — 48 (wl) 2 wl — 48 (wl) 2 wl + 9 6w k wlwl 
+100 (wl) 2 wl + 48 (wl)" wl — 64 wlw k w 7 k — 96 wlwlw 7 k — 48 (to ®.) 2 w k — 18 ( wl) 2 (to® + 4wl + 4w k 
+ 6 wl — 16 w 7 k + 24 w\) - 12toj((5 (wl) 2 - 3(10to)[ + 7w k + 4wl — 10 wl + 8 to®)to®. + 2(10 (wl) 2 

+2(8 w 5 k + 3 wl - 4w 7 k )wl + w 5 k (5w 5 k - 6 to£))))<| - 4c 3 (3toj[ - 3 wl + wl )(9 (wl) 2 - 3(3 wl - 5 wl 
+ 8 wl + 5 wl + 3 wl — 4 w 7 k )w\ + 5 (wl) 2 — 10 wlw k — 8wlwl + 6 wlwl — 6 w k w k + 6 wlwl + 3w k (—3w k 
+6 wl + 5 w k + 3 wl — 4 wl) + 4to|to|)c 2 + 3c 3 (3toj[ — 3 wl + to®) 2 (2tt4 — w k + to®. — wl — to®)] c\ 
+4cl(6wl - 3 wl + 2wl - 2wl) 2 [3wl + w 3 k - 4 (wl + wl - wl)}}, 

w 5 k+1 = + —^/^{16 [3 (wl) 3 + 16tof (wl) 2 + 14w 5 k (wl) 2 + 9wl (wl) 2 -46 (wl) 2 wl 

— 18 (wl) 2 wl — 80 wlw k w k — 24wlwlw k - 48 wlwlwl + 22 (wl) 3 + 48 wl (wl) 2 + 18 w k (wl) 2 
+36 (wl) 2 wl + 48w k w k w k + 9 (w k ) 2 (wl + 8 wl — 3to® — 16 w k ) — 8(3 (wl) 2 — (7w k + 3to®)tof 

+3 w k (wl + wl))w 7 k + 6wl(2 (wl) 2 + (12to)i + 4 wl + 3 wl — 20 w 7 k )w\ — 13 (wl)" - 9 (wl) 2 — 12 w k w k 
—4w k (8w k + 3 w k ) + 12 ( 2 wl + wl)w 7 k + 3 w k (w k + 3 w k + 6 wJ)) + 3 w k ( (wl)" + 2(4u^ + 3 w k + 8 w 7 k )wl 
+24wlw 5 k + 9( (w ^) 2 + (wl) 2 ) -4(7wl + 3wi)w 7 k )\c\ - 16[c 3 (l62 (w ^) 3 -9(15w^ - 26wl+40wl 
+24w\) (wl) 2 + 6(18 (wl) 2 - 2(26w k + 15wj>)w;j; + 9(w^ + w^)(3«;^ + w%) + 3w 2 k (-8w 3 k + 19w^ 
+ 9 w fc)) u; fc - 54 ( w fc ) 2 w l + 2(8 (u ^.) 3 - 3 ( llw k + Qw k) (wl) 2 + 3(w 5 k + w®)(llwfc + 3 wt)wl 
- 9w l( w l + w e k ) 2 ) - 3wl(l2 (w 3 k ) 2 - 4(13w k + 6 w%)w\ + 3 (w\ + wl)(llwl + 3^))) - 2c 2 (ll (wl) 3 
—64w k (wl) 2 — 114u^. (wl) 2 — 72w k (wl) 2 + 194 (wl)" wl + 72 (wl)" wl + 208 w k wlw k + 72 wlwlwl 
+234 wlwlwl — 80 (wl) 3 — 120 wl (wl) 2 — 72w k (wl) 2 — 144 (wl) 2 wl — 120 wlwlwl + 8(9 (u ;^) 2 

— (20 wl + 9 w k )wl + 9 w k (w k + w k ))w k + 18 (w ^) 2 (w^ + 2w 7 k ) + 9 (ui ^.) 2 (3 wl + 3u;[( — 40 wl — 36w^ 

—24u>® + 48 wl) + 3wl(3w 2 k (7wl + 20 wl + 15 wl + 9 w 6 k - 40 w 7 k ) + 2(7 (w ^) 2 - (52 wl + 67 w\ + 48w 6 k 
—60 w k )w k + 63 (wl)~ + 36 (wl)" + 84w^u>^ + 36 w k w k + 81u;fw® — 68 w k w k — 36w k w k )) 

+3wl(6 (' wl ) 2 + (24w^ + 27u^ + 21w® — 52 w 7 k )wl — 4(13 (w^) 2 + 17 wlw\ + 15w®w^ — 19 w k w k 
+9 (w®.) 2 + 3 wlwl — 9wlwl)))]cl — 32c2[c2(270 (wl) 3 — 9(9 wl — 68 w k + 96 wl + 140w^ + 114iy® 
—9Qw 7 k ) (' wl ) 2 - 6(9 (wl) 2 + 3(13 w\ - 36 wl - 46 wl — 42 wl + 44 w 7 k )wl — 63 (wl) 2 - 2(83 (w|) 2 
+8(15w® — 8 w 7 k )wl + 36wl(2w k + tuf) + 9wl(5w k — 4 wl)) + 6wl(20wl + 38w| + 29w® — 20^))^ 

— 18 (wl) (wl + 4^ + 4w k + 6 w® — 8 w k ) + 2(34 (w|) — (72 wl + 172^;^. + 117 w k — 72 w k ) 

x (wl)" + 2(119 K ) 2 + 88 wlwl + 156tt;®- 76 w k wl + 45 (wl) 2 + 36wlwl — 36 wlw k )wl 
-2w 5 k (47 (wl) 2 + 90 wlwl + 45 (w ^) 2 + 48 wl(w 5 k + ^)) + 72w 5 k (w 5 k + - 3w^(29 (wl) 2 

—4(22wl + 40 wl + 33wl — 28 w k )wl + 2(73 (w ^) 2 + 108 w k wl - 68 w k wl + 45 (w ®) 2 + 8wl(8wl 
+3 wl) — 36 wlwl))) — 603 ( 21*4 — 1 * 4 . + w k — wl — wl )(36 (wl) 2 — 3(9u>^ — lOit^ + 12io^ + 6 w k )wl 
+6 (wl) 2 + 6 (wl) 2 - 14wlwl - Qw k wl + 6w k wl + 3wl(—4wl + 7u| + 3io®))]ci + 3[64c2(6iUfc 
-3«4 + 2w| - 2w^)(2w^ -wl + wl-wl- wl)(3c 3 (-2wl + - wl + wjj + + 2c 2 (3w^ + 3w 3 k 

~4wl - 6 (wl + wl) + 4 ^))] }, 

= ^72 ~ ^J +72 { [ 216u, fc W ) 2 - 9 ( (^fe ) 2 + 4 ( 2w fe - 3w fc - + 6 w 2 k (wl + 2w r k ) - 4(2i*4 (! 

x (wfc + 6 w®) + 3( (wl) 2 + 3 w 6 k wl - 4w 7 k w\ - 6 wlwl)))w\ + 2( - 32 (to ^) 3 + 9(5w 3 k - 4w$) {wl) 2 
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-6 w\{wl - 10 wi)w\ - 27 (wlf + 27(wl - A(wl + wjj)) ( wif - 4 ( w 3 k f ( w 3 k + 6 w£) + 9(w 3 k - w%) 
x (w k + 16^4 + 13 w k )wl + 12(5 (wl) 2 - (11 wl + 15 w k )w k + 3 (w k + w%)(2w\ + 3 w k ))w 7 k ) - 9 wl 
x (( w lf + 4 0 1 + 3 wl + 4 w 7 k )wl + (wl) 2 + 24 wiwt + Aw 5 k (3w 6 k - 5 w 7 k ) - 9 w\(w% + 4 w 7 k ))\cf 
- 3 [ 3 c 3 (18 (wlf + 3(-30 wl + 15 wl + 16 (w\ - w 3 k )) (w 6 k ) 2 + 2(54 (wl) 2 + (-45 w 2 k + 66 w 3 k - 72 wl)w\ 

+9 (wl) 2 + (19wf - 21 wl)(w 3 k - wl) + wl(33wl - 30 w 3 k ))wl - (3wl + 2 wl - 3w 5 k )(6 (wl) 2 
+(-3 wl + 12 wl - 14 wl)w\ + A(w 3 k - w 5 k ) 2 + w 2 k (7w 5 k - 6 wl))) + c 2 (9(3w 2 k - 4(4 w\ + 8 w 5 k + 9w% 

— 16 wl)) (wl) 2 + 6(3wl(2wl + Awl + 5 w l + 3 wl — 24 w k ) — 2( (u>®)“ + (32 wl + 39 wl + 36 wl — 68 w 7 k ) 
xwl — 3(l6 (wl) 2 + 12 wlwl + 39wlwl + 15 (w k f + 28 w k w k ) + 4(19 wl + 27 w k )w k ))w k + 3w k (l5 

x (wl) 2 + 2(24 wl + 17 wl + 3wl — 60 w k )w k — 4(16 (■ w k ) 2 + (42w® — 34 w 7 k )wl + 14 wl(w\ + 3 w k ) 

+9w k (wf - 6 w 7 k ))) + 18 (wl) 2 (3 w 6 k + Aw r k ) + 4( (wl) 3 - (AOwl + A7w 5 k + A5wl - 6Aw 7 k ) {wlf 
+ (98 (wl) 2 + 88 wlwl + 225 w k w k + 99 (wl) 2 + 156 wlwl — A(3Aw k + 45 w k )w k )w k — 2(26 (wl) 3 
+90 wl (wl) 2 + 90 (w%) 2 w\ + 27 (wl) 3 + 6 wl(w 5 k + w & k )(Awl + 9w\)) + 36 (w 5 k + wl)(2w 5 k + 3wl)w 7 k ))] 
xcf + 12 c2 [c2(54 (wl) 3 — 9(3 wl — 14w| + 2Aw k + AOwl + 48 wl — AOwl) ( w lf + 6(19 (w k f — 2(28 wl 
+49u>jj + 51u>® — 38 w k )w k + wl(—3w\ + 30 wl + 45u>jj + 36 wl — 5Awl) + 2(41 {wlf + 30 wlwl 
+93 wlwl — AOwlwl + 45 (wl) 2 + A8w k w k — 54 w k wl))wl — 18 (w k f (wl + 2 (w k + w k — 2w k )) — 3w k 
x (6 (w 3 k f — A(13wl + 17 wl + 15 w k — 18 w k )wl + 65 (wl) 2 + 63 (wl) 2 + 56 wlwl + 96 wlwl + 144 wlwl 
—76 wlw 7 k — 108w k w k ) + 2(l6 (w^) 3 — (56 wl + lOlwjj + 105«4 — 68 w k ) (w k ) 2 + (157 (wl) 2 
+348 wlw k + 171 (wl) 2 + Awl(29w k + 42w®) - 20(7w^ + 9 wl)w k )w k - 3 (w k + w k )(Aw^(5w k 
+9 w%) + 3(8 (wl) 2 + 19w^4 ~ 8 w 7 k w 5 k + 9 (w^) 2 - I2w 6 k w 7 k )))) - 3c 3 (2w^ -w 2 k + w 3 k -wl- wl) 
x (18 (wl) 2 - 3(3«4 + 2(—5wl + 6u| + 9 w k ))wl + 3wl(—Aw\ + 5to| + 9 w k ) + 2(5 (w|) 2 
-(llw 5 k + 15w®)u>^ + 3(u>| + wl)(2wl + 3w®)))]ci - 4c^(6^ - 3 w 2 k + Aw 3 k - 4w^ - 6u|) [c 2 
x (72 (u4) 2 - 6(3 wl - 19 w 3 k + 2Awl + 37w^ + 33u4 - 24 w 7 k )wl - 9 (wl) 2 + 38 (w|) 2 + 92 (wl) 2 
+72 (w®) 2 — 30 wlwl + 72 wlwl — 72w\wl + 84w^w| — 130u>|w^ + 72wlwl + 72 w\w\ — 114u;^w;® 

+72WfcW® + 168 wlwl - 72(wl - w 3 k +w 5 k + w%)wl) - 27c 3 (-2«4 + w 2 k - w\ +w 5 k + wl) 2 ] }, 

w l +1 = { [ Hfe) 3 + (!6u4 + 18w| + 15w® + 58^4 - 192«;|) (w^) 2 + (96 (wl) 2 + 104^ u ,4 ( S11 ) 

+24u>®w^ — 304u4u4 + 45 (wl) 2 + 45 (wl) 2 + 168 (wl) 2 + ASwlwl — 372w k w k — 252 wlwl 
+6wl(Awl + 5u;jj + 9w® + 2Awl — 120w®.) + 432(w^ + w k )w k + 12wl(2wl + 3 wl + 24 w k ))w k 
+9 (wl) 2 (3 wl + 10 w 7 k - 72wl) + 3«4(48 (wl) 2 + 8(5 wl - 2Aw 7 k )wl + 19 (wl) 2 + 9 (w^) 2 + 112 
x (w 7 k ) 2 + 12w^w;® + 12 wlw 7 k — 204 w\w 7 k — 132w®w^ + 144(w^ + 2(wl + w k ))w k ) + 2( — 36 (wl) 3 
-78wl (wl) 2 + 158wl (wl) 2 - 63 (wl) 2 wl - 84 (w 7 k f w 5 k + 252 wlw 7 k w\ - 27 (wl) 3 - 84w^ (u>£) 2 
~ 72 ( w k) 2 ( w l + w k) + 90 ( wl) 2 wl +4 wl(- 23 (wl) 2 - 33wlw\ - 9 (wl) 2 + A2(w b k + w 6 k )w 7 k ) 

— 108 (wl + wl) 2 wl + 3wl(3 (w k ) 2 + Owlwl + 28w 7 k w\ — 28 (wl) 2 + 12tu®wJ + 4w^(2«;^ + 3w k 
+10 wl) - 72 (wl + w®)t«fc))]c? + [c 3 ((3wl + 2 w 3 k - 3 (w 5 k + w®))(27 (w^) 2 + 3(13w^ - 76 wl - 65w^ 

-33 wl)wl + 14 (wl) 2 - w 3 k (10Awi + 101 wl + 57u>^) + 12 ((wl + wl)(7w\ + 6w^ + 3w^) + w^(12^ 

+8 wl + 3^))) + 4(216 (wlf - 12(15^ - 19w^ + 21(^ + wl))wl + 27 (w^,) 2 + 59 (wl) 2 

+63 (wl) 2 + 63 (wlf — 102w^w| + 126w^w| — 126 w k wl + 126(u>^ — wl + wl)w k )wl) + c 3 (l9 
x (w^) 3 + (188w^ + 193?4 + 177w® - 832 w 7 k + 1200w^) (w|) 2 + 4( - 188 (wf) 2 + (-426^ - 282w£ 

+472 wl)w k — 261 (w®.) 2 — 189 (■u;®.) 2 — 252 (wl) 2 — 396u>;(«;® + 8lAw k wl + Q30w k w 7 k + 3w^(18wfc 
+10w^ + 3w® + A0w 7 k — 168to|) — 648(u>fc + wl)w k )w k + 9 (wl) 2 (7w k + 28wl + 37 w\ + 21 wl 
-232 wl + 480w®) + 6wl(l3 (;wl) 2 + (72w^ + 81 w 5 k + 57u>^ - 444w^ + 768 w 8 k )w 3 k - 2(112 (wlf 
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wl+i 


+2(115 w\ + 75 wt - 152 w 7 k )w\ + 133 ( w 5 k f + 2w\(9&w% - 239m;£) + 3(27io£ - 14ty^)(tyf - 4ty£)) 

+6 wl(6wl + 4m;®. + 3 m ;® + 22 w 7 k — 96 w k ) — 864(uj)) + w k )wl) + 4(209 ( w \) 3 + 495m;® ( m ;))) 2 
—600m;)) ( m ;)))“ + 405 (wl) 2 w\ + 252 ( m ;))) 2 w k — 1008w®w^w| + 135 ( m ;®) 3 + 252 m ;®. ( w |) 2 + 228 ( m ;|) 2 
x ( m ;)) + w k ) - 396 {wl) 2 w 7 k + 6wf(69 ( wl)~ + 110 'w k w k + 39 ( m ;®.) 2 — 84 (w k + wl)w 7 k ) + 324 (w k + m ;®) 2 
xw k + 9 (wl) 2 (w k + 12 wl) + 3m; 2 .(36 (wl)" + 4(13u;j) + 6 w% — 34 w 7 k )wl + 23 (wl) 2 + 18 (wl) 2 + 84 
x (w 7 k ) 2 + 30 w 5 k wl - 172 w\wl - 120 w & k w 7 k + 216(m;)) + wl)wl)))]cl + 2[(81 (wlf + 9(9m>)) 

+25 wl - 2(8 6wt + 107m;)) + 90m;)) - 218w 7 k + 264to|)) (w\) 2 + 3(71 (wl) 2 - 8(87 wf. + 110m;)) + 93w 6 k 
-194io£ + 204 w%)w\ + 2(288 (wl) 2 + 8(103u>)) + 78w 6 k - 80 w 7 k )w A k + 597 (w\f + 405 (w 6 k ) 2 + 336 (w 7 k ) 2 
+960 w k wl — 1300m;)) m;)) — 1068w®w^ + 86A(w k + w k )w k ) + 6w k (7w k + 6(4u^ + 5m;)) + 5 m ;®. — 22 w 7 k 


+40 wl)))wl + 65 (wl) 3 - 1348 (m;))) 3 - 864 (wl) 3 + 896m;)) (wl) 2 + 2178m;)) (w\f - 2208 w{ (wl) 2 
+1566m;)) (w k ) 2 — 1440w| (wl) 2 — 2808m;)) (m;®) 2 + 1008m;)) (w k ) 2 — 1008m;®. {wl) 2 — 1008m;® (wl)" 

—668 (wl)" wj. — 886 (wl) 2 wl — 960 (w k ) 2 w k + 2832 wlwlwl — 768 (w k ) 2 w k — 960 (w k ) 2 w k 
—3348 (wl) 2 wl + 2160 w k w k w k + 3552 wlwlwl — 3696 wlwlwl + 1348 (wl) 2 wl + 2904 (w k ) 2 w k 
+2088 (wl) 2 wl — 1952 wlwlwl — 4280 wlwlwl + 2016 wlwlwl — 3528 wlwlwl + 2016 wlwlwl 
+5040 wlwlwl + 36 (wl) 2 (2wl +wl + &(w 7 k — 4m;®.)) — 48(26 (wl) 2 — 54 (wl + wl)w\ + 27 (wl 
+w 6 k ) 2 )w 8 k + 3wl(ll (w 3 k ) 2 + 4(32m;)) + 37 w\ + 39w% - 126m;£ + 192 w%)w\ - 336 (w 7 k ) 2 - 4(64 
x (wl) 2 + 140 w k w k + 96 w k w k + 87 (wl) 2 + 63 (wl) 2 + 138w|w®) + 16(38«4 + 65 wl + 51 wl)w 7 k 
— 864(w^ + wl)wl))cl - 2c 3 (324 (wl) 3 - 9(21 wl — 67 wl + 96wf, + 145w^ + 105w® - 104u>£) (w^,) 2 
+9(41 (wfc)" — 2(50u^ + 82 wl + 61«;®)u^ + wl(—29wl + 84w;^ + 97 wl + 63 wl) + (w®. + wl)(108wl 
+127w 5 k + 81wg))u?J - 48(18^ - 20 wl + 21(«;® + w 6 k ))w 7 k wl + 74 (w®) 3 - 270 (wl) 3 - 162 
x (wl) 3 + 615«4 {wl) 2 - 252wl {wl) 2 + 405u^. (w^.) 2 - 252w| (u^.) 2 - 594w£ (wl) 2 - 232 (wl) 2 w^ 

-413 {wl) 2 wl + 492 wlwf.wl - 315 (wl) 2 wl — 702 (wl) 2 wl + 492 wlwlwl + 1020w^'u;®.w® 

—504u4u>jj'u;® — 3wl(30 (wl) 2 — (136w^ + 173wfc + 117 m;® — 152wl)wl + 6(m;^ + wl)(26wl + 26m;®. 

+15m;®. — 28 w 7 k )) — 18 (wl)" (8w^ + 6m;^ + 3 wl — 10 w 7 k ) + 4(61 (wl) 2 — 126(m;^ + wl)wl + 63(w®. 
+ w t) 2 ) w l ) c 2 + 3c^(2m>^ - wl + wl - wl - wl)(3w\ + 2 w 3 k - 3(w 5 k + wl))(9w 7 k - 6u^ + 4 m;^ - 3 
x (wl + w®))] ci + 4 c 2 [27cg(—2 m;^. + wl~w 3 k +w 5 k + wl) 3 + 9c 2 c 3 (21m;^ + 15m;^ - 4(7u;^ + 9 m^ 

+6u;® — 7w k ))(—2w k + wl — wl + w k + m;®.) 2 + c 2 ( — 270 {wl) 3 + 9(9m;^ — 69u^ + 128 m;| + 189m;®. 

+165u>® — 224 m;J + 192m;®.) (wl) 2 + 3(9 (wl) 2 + 6(11m;| — 36m;^ — 50it;®. — 48 m;® + 84m;J — 96m;®.) 
xw 2 k + 2( - 74 {wl) 2 + (236m4 + 361m^ + 309m;^ - 380u>J + 288 wf)wl - 2(56 (m>^) 2 + 202^.^ 

+162u;®m4 + 175 (m;®.) 2 + 126 (w 6 k ) 2 + 56 (w 7 k ) 2 + 294m;^m;^ - 2(56m4 + 137m;| + 117 w%)w 7 k ) 

—288(m;®. + wl)wl))wl - 101 (w 3 k ) 3 + 656 ( W ®.) 3 + 432 (m;®) 3 - 336m;® (wlf - 1224m;® (m;®) 2 
+920m4 (wl)" — 900u>^ (wl)" + 648m;)) (m;®) 2 + 1440m;| (m;|) 2 — 336m;®. (wl) 2 + 336m;®. (wl) 2 
+336u>® (wl)" + 416 (wl) 2 wl + 669 (wl) 2 wl + 336 (wl) 2 w k — 1336m;))m;))m;®. + 573 (wl) 2 w k 
+336 (wl) 2 wl + 1680 {wl) 2 wl — 1080 wlw k wl — 2064m;®.w®.m;® + lh8Awlwlwl — 632 (wl) 2 wl 
-1136 (wl) 2 w 7 k - 864 (wl) 2 wl + 672 wlwlwl + 1768 wlw 3 k w 7 k - 672ui^m;®.m;^ + 1512m;®.m;®.m;£ 

— 672wlwlwl — 2016m;®.m;®.m;| + 432(—m;)) + m;^ + m;®) 2 m;® + 9 (m; 2 .) 2 (m;®. + 4(ui^ + 2 m;®. + 3m;®. — 7 wl 
+12m;®.)) + 3m; 2 .(31 (m;^) 2 - 4(38m;^ + 53 w\ + 48m;£ - 74m;^ + 72 w%)w\ + 4(28 (wlf + (80m;)) + 60m| 
-56w 7 k )wi + 61 (wl) 2 + 45 (wl) 2 + 28 (m;£) 2 + 102m;)>)) - 116m;^m;^ - 96w e k w 7 k + 72(m;)). + m;^)m;|)))] }, 

<^72 + 27 p 9/2 {[ 112 ( w fc) 3 + 36w fc ( w fc) 2 + !32 m;^ (wl) 2 - 36u>® (m 4) 2 - 336u>^ (wl) 2 (S12) 


+78 (wl) 2 wl + 54 (wl) 2 wl + 336 (wl)" wl + 108m4u>®m;! 


— 36 wlwlwl + 2\Qwlw 7 k wl + 72 wlwlwl 



6 


—624 wlw 7 k wl — 288 w%w\.w\ + 27 (wff + 27 (wf) 3 — 112 (wf 3 + 81u>| (wff — 468wj[ (wf 2 
+108w^ (wf 2 + 528u>jj (wf 2 + 360w® (wf 2 + 27 (wf 2 w® + 54^^^® — 348 (w k f w 7 k — 216 (w®) 2 
xu>l + 216 w\w\w 7 k + 3 6w k w k w k + 108 w\w k w 7 k — 432w k w k w k + 108(3 (wf 2 — 2(8 w k + 7w k + 3 w k 
-8 wfw\ + (w k + wf(8wl + 7 w k + 3il ;®. - 8 wf + 4wl(2w k + wl - 2wf)wl + 6 (wf 2 (wl + 2.w 7 k 
+24 wl) + 6w^(2 (wf~ + 8 w k wl + 6 w k w k + 28wlwl + 2 (wl) 2 — 46 ( w k ) 2 + 3 w k w k + &w],w 7 k 
+24 wfw 7 k + 18 wlwl + 72(wl — 2 w k — 2wl — wl + 2wfwf] cf + 3c 3 [9(10w^ + 7 wl + 3u>® — 24 w k 
+48 w k ) (wl) 2 + 6( - 48 (wl) 2 + 12 (wl + 8w^ + 9 w\ + 5 wl)w 7 k - 3(16 (wl) 2 + 26 wlwl + 11 (u>jj)“ 

+3 (wl) 2 + 10(«4 + wfwf) — 72(wj? + wjj + w k )w k + 2w^(10u4 + 7w k + 3w® — 22w| + 36 wf)wl 
+4 (w®.) 2 (10«4 + + 3u>® — 20 wl + 27u>®) + 12 w k ( — 12 (wl) 2 — 2(1 lw| + 5w®. — I2w 7 k )wl 

— 10 (wl) 2 — 3 (wj®) 2 — 12 (wf 2 — 10 w\w% + 4 w k w k + 30 wlwl + 18 w k w k — 18 (w k + w\ + wf)wf 
+9(12 w 8 k (w 2 k ) 2 + 4(4 (wfc) 2 + 4w 5 k wl - 8w£u>£ + (wl) 2 + 4 (w|) 2 - 6 w%w 7 k - 2\wfw 7 k + Q(w\ 

+wf)wl)wl + (wjj + W®)(16 (wl) 2 + 26w^w^ + 10 w%wl - 32 w 7 k w\ + 11 (■ w 5 k ) 2 + 3 (w 6 k ) 2 + 16 (w 7 k ) 2 
+10 wlwl - 32wlw 7 k - 16 wlw 7 k + 12 (wl + w®)w®))]cf - 9c§(2«4 - wl + wl~w 5 k - wl) [9 (wl) 2 
+6(2 w\ — 8 wl — 7wl — 3 wl + 8 w 7 k )wl + 4 (wf~ + 3 (w k + wf)(8w k + 7w\ + 3 wl — 8wl) 

+12wl(2wl + wl — 2 wl) — 12wl(2wl + 2wl + u>® — 2toJ)] c\ + 27c\(—2w\ + io(! — tu 3 + + w®) 3 

+2c 3 [27 (wl) 3 + 27(3 w 2 k + 7w\ - 2(22 w 4 k + 32 W ® + 33 W ® - 70^1 + 96^)) K) 2 + 9(21 («; 3 ) 2 
—4(46 wl + 6 6wl + 63w® — 130wJ + 168 w k )w k + 16(13 (wl) 2 + 37 wlwl + 30 wlwl — 50 wlwl 

+27 («4) 2 + 18 (wl) 2 + 37 (w k ) 2 + 45'u;®.w®. — 76 wlwl — 57 wlwl + 6(8 wl + 13w| + 9u>® — 8 wl)w k ) 

+6 wl(wl + 4(wl + 2w^ + 3w® — 7w r k + 12 w k )))wl + 51 (w^) 3 — 448 (wl) 3 — 1824 (wl) 3 — 864 (u>®) 3 
+448 (w 7 k ) 3 + 1104it> 3 . (w^) 2 + 2376«4 (wl) 2 - 3408m)^ (wlf + 1512 w\ (wlf - 2160^ (wff 
-3456 wl (w 6 k f + 2832w\ (w 7 k f - 1344wl (w k f - 3840w^ (wlf - 3168w®. (wl) 2 — 564 (u’lf wl 
-792 (w 3 k f w 5 k - 2112 (wlf wl + 3216 wlwiw\ - 702 (wlf wl - 1440 (wlf w 6 k - 4320 (w 5 k f wl 
+2592 w\w\wl + 3888wlwlwl — 5184 wlw k w k + 1428 (wlf wl + 1344 (wlf w 7 k + 5568 (w®.) 2 
+3456 (wlf wl — 3936iulwlwl — 6240 wlwlw 7 k + 5952 wlwlw 7 k — 4752 wlwlw 7 k + 4608w^w®w| 
+8640 w^m;®w’( — 864('ui 3 . — u;®. — wl)(2w\ — 4w^ — 5w^ — 3w® + 4u^)'u; 3 + 9w^( (w^)“ + 8(3u>[[ 

+5 wl + 6w® — 15 wl + 24w®)'ui 3 . — 240 (w k f — 8(6 (w k f + 18(w| + wl)wl + 13 (u>^) 2 + 9 (wlf + 24wjj 
xw[) + 96(3 wl + 4wfc + 3wl)w k — 96(4 wl + 5w®. + 3w k — 4w k )w k )] + 12c| [ci ((wlf + (34u;^ + 36w^ 
+33w® — 170 wl + 336 wl) (w k f + ( — 160 (w k f — 364 w k w k — 288 wlwl + 712 wlwl — 218 (w|) 2 

— 135 (u>®)“ — 564 (wfc)” — 348 w k w k + 944«;^wJ + 684 wlwl + 6wl(2wl + wl + 14■uij — 48w 3 ) — 432 
x + 3 wl + 2u>® — 2 wl)w k )wl + 112 (w^) + 232 (u>|) + 108 (wl) — 112 (w[0 + 24w^. (wl) 

+42wl (wlf + 519w k (wlf + 27w^ (w®.) 2 + 306^^. (wlf + 432w^ (wlf + 396 wl (wff + 336 wl (wff 
+816w| (wf 2 + 648w® (wf 2 + 396 (wf)~ wl + 72 wfvlwl + 228 (wlf wl + 522 (wl)~ wl + 108u>^ 
xw k w k + 90 w k wlwl + 708 w k wlwl — 336 (wf"w 7 k — 972 (wlf wl — 576 (w®)“ wl — 408 w k w k w 7 k — 420 
xw k wlwl — 1200 w k wlwl — 288 w k wlw 7 k — 864 w k wlw 7 k — 1440^;^;®^^ + 108((w^ + wf(8wl + 9w® 
+5w® — 8w^) + 4wf2wl + 2w®. + w® — 2 wf)w k + 9 (wj.) 2 (iy^ + 4w^ + 7w®. + 6u>® - 46u^ + 108w®) 
+3wl(2 (wlf + 2(12 wl + 15 wl + 15w^ - 88w 7 k + 192^)u>^ - 84 (wff - 113 (w^) 2 - 72 (wf) 2 

—348 (w 7 k f — 188 wlwl — 144wlwl — 174 w\w\ + 424w^w^ + 536 wlwl + 372 w k w 7 k + 6wf2wl + 

+6^1 - 24w^) - 72(8w k + llw| + 7u>^. - 8 wfwf)) - 9c 3 (2wl - w\ + - wf)( 9 (w^) 2 + 4(3 wl 

—7 wl — 9w k — 6w® + 10wJ — 6 w k )wl + 4 (wf 2 + 16 (w^) 2 + 28 (wf) 2 + 12 (w®) 2 + 16 (w 7 k f + 2 w^wl 
+3 w k wl + 42 wlwl + 3wlwl + 26wfwl + 36 w^w® — 8 wlwl — 32w k w 7 k — 48wlw 7 k — 32wlw 7 k + 12(u>^ 
+w^ + wf)w\ — wl(18wl + 23u>jj + 15w®. — 24w^ + 12u> 3 ))] + 3c2 [(9(3w® + 22 wl — 144w®) (wf 2 
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+6(20 ( wlf + 40 w k w k + — 248 w\.w\ + 23 (m;®.) 2 + 18 (m;®.) 2 + 252 (wl)" + 24m;®.u;®. + 6u;j: 

xw 7 k — 256wlwl — 156 w k w k + 3wl(4wl + 2w\ + m;® + 14u;J — 88 w%) + 72(w k + 8wl + 9w k + 5 u;® 

-8 w 7 k )wl)wl + (w 3 k y (24+ 12m;| + 3m;® + 94m;£ — 480m;®) + 2w\(b2 (wlf + (80 w\ + 60tu® 

—440u;J + 864 w k )wl + 43 (m;|) 2 + 8 w k (9w k — 59wJ + 135 m;®) + 6(6 (it;®.) 2 — 54 (m;J — 2 w k )w k + w k (w k 
+70wl) + 24(w 2 k - 6wl)wl)) + 2(112 (m;£) 3 - 168(2^ + 4m>£ + 3 wl) (w 7 k f + 24(14 (m;£) 2 + 38 \w\w\ 
+24w 6 k wi + 26 (wl) 2 + 15 (wlf + 36 w b k wt)w 7 k - 2(56 (w^) 3 + 12(11^ + 4 m;£) (w^)“ + 6(21 (m;|) 2 
+22m;®.m;® + 12 (wt) 2 )wi + 49 (wlf + 27 (m;®) 3 + 99 w 5 k (w 6 k f + 99 (w 5 k f w 6 k ) - 432(m;® + m£)(2u,£ 

+2wl + wl — 2w k )w k + 3wl(l2 {wlf + 12 w k w k + 56 w k w k + 3 (w k f — 84 ( w k f + 40 w\w 7 k + 24io®m;| 
—72(4 w k + 3 wl + to® — 4 w k )w k )))c 2 + 4cs(27 (wlf + 9(6 w k — 32 w k — 29 wl — 15m;® + 52 wl — 48 w k ) 
x (wlf + 3(12 (wlf - (114t0fc + 107w| + 57u;j? - 176 m;£ + 144w^)u;f + 3(32 (w^) 2 + Uw 2 k w£ 

+82 w k w k + 50m;® m;£ — 64 w k w k + 48 (w k f + 18 (w k f + 32 (w k f + 9wlwl + 3u^u>® + 58 m;)? m;® 

—32 w^wj. — 104 wlwj. — 72wlw\ + 48 (m;^ + m;)? + wl)wl))wl + 8 (wl) 3 — 2 (wlf (50wl + 49m)J? + 27u>® 

—74 wl + 54m;®) + 3 m;®. (48 (wlf + 28 w k w k + 134m ;®.wl + 86m;® tu^ — 9 6w k w k + 83 (m^)“ + 33 (m;®) 2 
+48 (m;^) 2 + 18m;^m;| + 6m;^m;®. + 104m;®. m;® — 56 wlw 7 k — 168m;®.m;£ — 120 w k w k + 72(m^ + wl + wl)w k ) 
+9(4(m > 7 k - 3w^) (wlf - (16 (m;|) 2 + 30wlwl + 14wgw;£ - 32 m>|m>)‘ + 13 (w 5 k f + 3 (m;£) 2 + 16 (w 7 k f 
+12m;®.m;® — 40m;|u;£ — 24 w k w 7 k + 24(u;)j + w k )wl)wl — 2(wl + m;®)(8 (m;^) 2 + 17m;^m4 + 9 wfwl 
— 16 w k wl + 9 (m;|) 2 + 3 (m;®.) 2 + 8 (w k f + 10m;|m;® — 2Qwlw 7 k — 12m;^m;^ + 6 (w k + m;®.)m;®)))ci + 36 
xc 2 (—2wl + w\ — w k + m;®. + wl) 2 (3wl + 2wl — 4wl — 5m;)? — 3m;®. + 4 w k )] }. 

It is straightforward to show that in the zero-held limit s k = r k = r k = wl = w k = ■ ■ ■ = wl = 0, and Eqs. (S2), 
(S5) reduce to the RG equations for the hierarchical Edwards-Anderson model [2]. 


S2. RENORMALIZATION-GROUP FLOW IN TERMS OF PHYSICAL PARAMETERS OF THE MODEL 


Here we will discuss how to express the RG flow of the coefficients 


K>, Kl 


(S13) 


in terms of the physical parameters of the model, i.e. the interaction-decay exponent <r, the inverse temperature f3 
and the magnetic-field strength = E [h 2 ]. While the RG FPs can be completely determined by imposing that the 
FP is invariant with respect to the RG transformation, to determine the RG coefficients for any finite k an explicit 
solution of the RG recurrence Eqs. (S1)-(S12) is required. Given that the solution of these equations is particularly 
involved, we will not provide the explicit expression of s k , {r^}, {m;£} as functions of a, /?, <; k here, but we will rather 
lay out the procedure to determine it. First, we will relate the coefficients (S13) to the averages of moments of the 
overlap Q a b —see for example the RHS of Eq. (19). Second, we will write the average moments of Q a b for k = 0 in 
terms of a, (3 and <+: as a result, we will obtain an implicit expression of sq, {tq}, {m>q} as functions of the physical 
parameters of the model. The dependence of the RG coefficients (S13) on a, 13 and g, can then be obtained implicitly 
from So, {?’g}, {^’o! by iterating the RG Eqs. (S1)-(S12). 

In order to illustrate this procedure, let us introduce the indexes A, B,C ,..., which run over replica pairs 
(1, 2), (1, 3),..., (n — 1, n). In addition, let us introduce the coupling constants 


p _ < d 3 flQ] 
9abc ~ 6 dQ A dQ B dQ c 


(S14) 


By using Eqs. (7), (S14), the replicated partition function can be rewritten as 


2 s k ^ ( Qa 


C . 

T 5Z Qa*abQb 

AB 


1 
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E&bcQaQbQc 

p =1 


Zk[Q\ = exp 


(S15) 



Equation (S15) has the structure of a field theory with a quadratic term and a non-quadratic correction: this form 
allows for extracting easily the overlap moments by means of diagrammatic methods [3]. The result is 


Qa = 

QaQb = 


4 s k Jf 


-pE E sIcd^a'b 


p =1 B,C,D 


— (2 s k Jif) 2 + ^CD 


0(5 2 ), 


(516) 

(517) 


QaQbQc = — 


(^) + ^xh + ^^iz&DE'PXc^h + oig 2 ), 

^ ' p —1 

(^) 3 - {&AB + ^AC + ^c) - (§) 3 &Zb&Bh&5b + o(g 2 ), (S18) 


where ^ = £" <6=1 ^ 12^6 

We will now write the LHSs of Eqs. (S16)-(S18) in terms of Boltzmann averages of overlap products (qu), {$ 2 ), 
..., where in what follows the average () of a function of multiple replicas will denote the average over all replicas, 
each replica having an independent Boltzmann measure with Hamiltonian (1). Let us start with Ecp (S16): we have 


E[(<?i2}] = lim f[dQ\Z k [Q\Q 12 (S19) 

n^Oj 

= lim C -(^+i)f f [dQ]Z k [Q]Q 12 , 
n^O J 

= C~ k ' 2 lim 0^2; 

n—t 0 

where in the first line of Eq. (S19) we used Ecp (2), in the second line we used Eq. (5), and in the third line we 
utilized the normalization identity lim rw o / [dQ] Z k [Q\ = 1 to write the integral in the second line as an average with 


respect to Z k [Q\, which we denote by . Proceeding along the lines of Eq. (S19), we obtain 

E[(q A q B )} = C~ k lim QaQb, (S20) 

n —>-0 

E [(q A q B qc)] = C~ 3k/2 lim QaQbQc ■ (S21) 

n —>-0 

For k = 0, we can write the RHS of Eqs. (S19)-(S21) as functions to the physical parameters /3, ft: 

E[(< 7 i 2 )] = E[(qi 25 i 3 )] = E[(gJ 2 )] = E[(^ 2 gi 3 )] = E[(^ 2 g 34 )] = E[(q 12 q 13 q 2 i)} = E[tanh(/3/i) 2 ], (S22) 

EK 912 )] = E[(gi2<?23 ( ?3i)] = b (S23) 

^[( 912934 )] = E[(gi 2 (?i 3945 )] = E[tanh(/3/i) 4 ], (S24) 

E[(gi29i39i4)] = 0, (S25) 

E[(gi 2 g 34956 )] = E[tanh(/3/i) 6 ], (S26) 


where the expectation values of powers of tanh(/i) in Eqs. (S22), (S24), (S26) depend on /3, ft only through the 
product /3ft. 

By using Eqs. (S16)-(S18), (S19)-(S21), (S22)-(S26) we obtain a set of twelve equations which relate the RG 
coefficients (S13) for k = 0 to a and /3ft: 

s 0 = s 0 (cr,^ft), rg = rg(cr,/3ft), wg = wg(cr,/3ft), (S27) 

where in the LHSs of Eqs. (S27) we imply the n —i 0 limit. The RG coefficients for any value of k can then be related 
implicitly to a and /3ft by iterating of the recursive Eqs. (S1)-(S12) with initial condition (S27). 
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